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Abstract 

A threshold singularity in a heavier meson channel, e.g. D*D*, generally reflects 
in the amplitudes of e~^e~ annihilation into lighter meson pairs, such as D*D, DD, 
and DgDs, due to the channel coupling. The behavior of the cross section of the e+e~ 
annihilation into each channel in this situation is considered using general properties 
of the amplitudes near the new threshold. It is argued that the effects of the channel 
coupling can be quite significant in the observable cross section, in particular near the 
D*D* threshold. 



It is quite well known[l^ that the cross section of e^e^ annihilation in the region just 
above the threshold of open charm production around 4.0 GeV displays an intricate behavior 
which is yet to be studied in detail. This behavior is contributed by the successive onset 
of specific channels with the D mesons: DD, D*D, DgDs, etc., by the strong dynamics in 
each of these channels and by the coupling between them. Thus a detailed study of this 
region, as well as of the similar region at the B^*^ meson threshold, may provide a wealth of 
information about the strong dynamics of systems with heavy and light quarks. It has been 
recognized long ago [21 01 IH IHj that such dynamics can be quite interesting due to the fact 
that the strength and the range of the strong forces between the heavy mesons is fixed by 
the interaction of the light quarks and gluons, while the mass of the interacting constituents 
is set by the heavy quark mass. In particular this behavior suggests existence of bound 
or resonant states of "molecular" type [2] at least in some meson-antimeson channels. A 
recent discovery 6^ of the resonance X(3872) and its subsequent studies[ZllHlinilini[II] niost 
strongly support the interpretation[T21 [1311111121 of this resonance as a "molecular" 5" wave 
state in the D^D*^ + D^D*^ channel with quantum numbers J^'" = 1^^, which however is 
inaccessible in the e~^e~ annihilation. 

It is quite apparent phenomenologically that the P wave states of the meson pairs pro- 
duced in the e'^e" annihilation have a nontrivial strong dynamics of their own. In particular 
both the DD and BB thresholds start with the all well known iIj{3770) and T(4S') resonances, 
and there is a long known conspicuous enhancement of the e~^e~ cross section near the D*D* 
threshold just above 4.00 GeV, which has been discussed long ago||2lj as a "molecular" D*D* 
resonance. Thus existence of a resonance near the threshold of a specific heavy meson pair 
looks more like a rule than an exception, and it is quite possible that similar singularities 
will be observed at other thresholds, such as D*D, B*B, etc., once a more detailed study of 
the e~^e~ annihilation at those thresholds is done|16j. 

Admittedly, the current theoretical understanding does not provide much of an insight 
into either the properties, or even the very fact of existence of the threshold resonances. 
Furthermore, some of the usual guidance from the heavy quark limit is also lost in this 
situation. An example of such loss is the inapplicability of the conservation of the heavy 
quark spin at the energy scale relevant to the discussed resonances. Indeed, the part of the 
Hamiltonian, which depends on the chromomagnetic interaction of the heavy quark spin, 
is normally suppressed by the inverse of the heavy quark mass tuq and is of the order of 
Ag^^/mg, which is also the scale for the mass splitting between the vector and pseudoscalar 
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heavy mesons. However at energies in the region of the thresholds of combinations of the 
vector and pseudoscalar mesons, the relevant overall energy scale is set by same mass split- 
ting. Thus the spin-dependent Hamiltonian cannot be considered as a small perturbation in 
this region. This behavior in particular invalidates in this energy region any predictions of 
the relative yield of meson pairs with mesons of different spin, based on the heavy quark spin 
symmetry, as was demonstrated 3, in the discussion of a threshold resonance in the D*D* 
system within a simple spin correlation model, which model is equivalent in this channel to 
a more general scheme used in a coupled-channel model of charmonium|17| ITSj . It should 
be also mentioned that in the latter approach the breaking of the spin-counting rule and 
the structures in the cross section of the e^e~ annihilation are describedllT^ in terms of the 
coupling between the charmed meson pairs and the charmonium dynamics rather than in 
terms of resonances in the charmed meson pairs, which approach is alternative to the one 
considered here. 

In this situation it appears to be worthwhile to discuss the behavior of the production 
amplitudes across the specific meson-pair thresholds in a model independent way, using the 
general quantum mechanical properties of unitarity and analyticity. Although such approach 
inevitably involves phenomenological parameters, which currently can only be determined 
from data, this looks like a necessary step toward understanding the strong dynamics of 
the heavy hadrons. Furthermore, in view of a significant and rapid variation of the e'^e" 
cross section in the threshold region, a proper parametrization of the cross section is of a 
practical importance for unfolding the radiative corrections in a concrete detailed scan of 
this energy region, such as being undertaken in the CLEO-c experiment |19j . In this paper 
the behavior of the amplitudes for production of a lighter and a heavier meson pairs is 
considered immediately below and above the threshold of the heavier pair. For definiteness, 
the discussion will refer to the threshold of the 'heavier' D*D* pair production, where the 
kinematically available 'lighter' channels are DD, D*D, and DgDs, and where the situation 
is likely to be interesting due to apparent existence of a threshold resonance, although the 
discussed properties may also be relevant for other similar thresholds. It will be shown that 
a reflection, due to the channel coupling, of the D*D* threshold singularity in the lighter 
channels should give rise to a somewhat nontrivial behavior of the cross section in those 
channels. It should be mentioned that the reflection of the threshold onset in one channel on 
another coupled channel is well known (see e.g. in the textbook [2111, Sect. 147). However 
the specifics of the pure P wave dynamics in the meson production amplitudes in the e^e~ 
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annihilation and of the parameters pertinent to the D^*^ mesons (and possibly to B^*^ as 
well) bring in interesting features, which justify a separate study of this phenomenon in the 
discussed context. 

In what follows the variation of the production amplitudes is considered in a rather 
narrow energy range of energy immediately below and above the D*D* threshold, on the 
order of about 30 MeV on each side. In this interval the D* mesons are nonrelativistic 
since their momentum k (real above the threshold, and imaginary below the threshold) does 
not exceed approximately 250 MeV. Also such momenta are reasonably small in comparison 
with the inverse of the distance scale of the essential strong interaction, and the expansion 
of the dependence of the considered amplitude on the D* momentum can be limited to 
the first terms mandated by the general quantum-mechanical properties of the amplitudes. 
The lighter mesons, on the other hand, are already quite fast at the energy of the D*D* 
threshold, e.g. the momentum p of each of the mesons in the D*D pair is about 500 MeV, 
so that it cannot be considered as small. However the variation of this momentum over 
the considered narrow energy range is relatively weak, so that the proper dependence of the 
amplitudes on p can be neglected in the first approximation, and we can rather concentrate 
on the dependence on the momentum k whose variation is quite significant. 

There is a minor complication of the behavior in the threshold region due to the isotopic 
mass difference between the charged and the neutral D* mesons, which amounts to existence 
in fact of two separate thresholds split by about 7 MeV, and this splitting will be accounted 
for in the discussion here. Another potential complication arises from the Coulomb interac- 
tion between the charged D* mesons, described by the parameter a/v, with v being the cm. 
velocity of the produced mesons. The effect of this interaction can be quite substantial at 
energies of about 5 MeV on each side of the threshold for the charged mesons. However, the 
Coulomb effects in the presence of strong interaction, especially in a situation, where there 
is a nearby resonance, have a peculiarity of their own |2H 1^ and merit a separate study. In 
the present discussion these effects are ignored, so that the resulting formulas do not contain 
finer details, which may be present very close to the charged meson threshold. 

We first start with considering the effects of the channel coupling ignoring the isotopic 
mass difference between the D* mesons, and will later modify the resulting formulas by 
including the mass splitting effects. It should be also mentioned that, although for the 
mesons produced in the e~^e~ annihilation there is only one relevant partial wave, namely 



3 



the P wave^, there are still several channels to be considered in the most general treatment. 
Namely, the 'light' channels are D*D, DD, and DgDs, and for the opening at the considered 
energy channel D*D* there are still two possible states differing by the total spin of the vector 
meson pair: 5 = and 5 = 2. The treatment of the coupling between the channels greatly 
simplifies in the presence of a resonance. In this case each of the production amplitudes 
Al in the light channels can be approximated as a (complex) constant ai, the non-resonant 
background, and the contribution of the resonance, through which the coupling to the 'heavy' 
channel takes place. As to the 'heavy' channel itself, the non-resonant production is small 
near the threshold due to the P wave kinematical suppression, and the amplitude can be 
approximated by that dominated by the resonance. The resonance couples to a specific linear 
combination of the two possible spin states, and the production of only this combination is 
relevant in the threshold region^, and the corresponding production amplitude for the D*D* 
pair is denoted here as Ah- 

Let W denote the total cm. energy E in the e~^e~ annihilation relative to the threshold 
of production of the two vector mesons, each with the mass M: W = E — 2M, and let the 
'nominal' position of the resonance be at the complex energy corresponding to W = Wq — 
iro/2. According to the general quantum-mechanical consideration, based on the unitarity 
and analyticity of the production and scattering amplitudes (see e.g. in the textbook pUj . 
Sects. 133 and 145) at energy above the 'heavy' threshold, the resonant production amplitude 
A2 has the form 

A = biA (u 

and the production amplitudes for each of the light channels are parametrized as 



where k = yMW is the cm. momentum of the D*D* pair. The coefficients a^, b^, and are 
complex parameters, which are also taken to be constant in the present consideration. The 
magnitudes and the complex phases of these coefficients result from the rescattering among 
the 'light' channels and the corresponding absorptive parts associated with these channels. 



^Generally an F wave production amplitude is also possible for a pair of vector meson. However in the 
threshold region this higher partial wave can obviously be neglected. 

Clearly, this reasoning rather naturally assumes that there is only one ' resonance in or close 

to the considered energy region. The spin structure of the state of the D* mesons at the resonance can be 
determined experimentally from angular correlations 
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It is worth noting that there is no reason to expect the phase of the coefficient bn to vanish 
at A; = 0, due to the existence of 'hght' inelastic channels strongly coupled to the 'heavy' 
meson pair. The normalization convention used here for the amplitudes corresponds to the 
cross section a{e~^e~ — > H) = {AhI"^ k for the 'heavy' channel and a{e~^e~ —>■ L) = \Ai\'^pi 
for each of the 'light' channels with being the cm. momentum in the corresponding 
'light' channel at the 'heavy' threshold, which momentum is taken as a constant across the 
considered energy range. 

The expressions and (j2)) can be readily found from the standard summation of the 
'blobs' in the graphs of the type shown in Fig.l, where the ultraviolet part of the loop 
is absorbed into the definition of the overall normalization and of the position Wq of the 
resonance. However the general treatment of the threshold behavior of the amplitudes |2Uj 
guarantees that these expressions are quite independent of assumptions about the form 
factors used in the graphs of Fig.l. 



L 01 H L 01 H L 01 H 



Figure 1: The graphs for the self energy of the resonance X near the threshold for the 'heavy' 
channel (H). The absorptive and dispersive parts of the self energy arise from loops with 
the 'light' (L) and heavy meson pairs. The ffiled circle stands for the electromagnetic vertex. 



Furthermore, the 'nominal' width Fq is due to the coupling to the lighter channels, and the 
term gjj in the denominator in equations and represents the additional contribution 
to the resonance width associated with the decay into the 'heavy' meson pair, and g is the 
coupling constant with the dimension of length. This extra contribution to the width can 
be written in terms of the excitation energy W as k^ = W ^Jw/w, where 

w = g~^M-^ (3) 

is a parameter with dimension of energy. The critical point of the present consideration is 
that the amplitudes Al as functions of the energy W are analytical functions of W (with a 
cut at positive W) and the expression can be analytically continued to negative values of 
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W. Using the notation A = —W at W < 0, the amphtudes Al below the 'heavy' threshold 
take the form 



The term in the denominator in the latter formula, proportional to A^/^, is clearly the 
analytical continuation of the g^k"^ term in eq.Q, and the proper sign is as indicated in 
eq.(0])^. This term becomes more essential than the one linear in A, starting from A ^ w, 
and gives rise to the effects discussed in the present paper. In particular, the expression 



A = Am = 16w/9. Thus in the case where the resonance is below the D*D* threshold, i.e. 

at negative Wq = —Aq, the real part of the resonance denominator has either two zeros if 

Ao < 16w/27, or no zero at all otherwise. In the latter case the position of the minimum of 

the absolute value of the denominator, i.e. of the maximal effect on the cross section, is at 

A = Am and does not depend on Aq at all, but rather is determined by the parameter w. 

The energy corresponding to A = A^ in fact also determines the position of the maximum 

of the resonance amplitude in the case where the 'nominal' position of the resonance is 

above the threshold, i.e. at positive Wq. This is due to the fact that at sufficiently small 

w the damping of the 'light' channel amplitude due to the absorptive term gj^ in the 

denominator in eq.© suppresses the 'light' channel above the 'heavy' threshold. 

Clearly, the significance of the discussed threshold effect critically depends on the value 

of the parameter w for the near-threshold resonance. At least for two known resonances, 

ipiSTTO) and T(4S'), at the thresholds of pairs of pseudoscalar heavy mesons the value of 

w can be found from the data: w ~ 60 MeV for ?/^(3770), and w ^ 20 MeV for T{AS). 

There is every reason to expect that for the possible resonance near the D*D* threshold the 

coupling of the resonance to the vector meson pair should be stronger than for the case of 

pseudoscalar mesons, and the parameter w should be significantly smaller than for iIj{3770). 

Indeed, according to the datapP, the observed peak at the D*D* threshold^ is mostly due to 

a very strong production of the D*D* pairs, and it has a width of 50 ± 10 MeV, in spite of 

a substantially smaller momentum k of the vector mesons as compared to the momentum 

■^This sign is due to the P wave motion, and is opposite to the one that would arise in a similar term for 
an S wave (ex — A-'^/^). 

''This peak is referred to in the Tables |Ij as V'(4040) with the maximum of the total e^e~ annihilation 
cross section at 4040 ± 10 MeV. However due to the threshold effect the maximum of the total cross section 
does not correspond to the position of the resonance, which position is in fact at a lower energy. 



Al = aL + 




(4) 



|Ay^— A cannot be arbitrarily negative and reaches its minimum, equal to —lQw/2'1 at 
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of the D mesons at the peak of '?/^(3770). Thus the couphng g for the vector mesons to 
their threshold resonance should be significantly stronger than of the pseudoscalar ones to 
?/'(3770). The parameter w is proportional to g~'^ (cf. eq.®)- Therefore, it is quite likely 
that for the D*D* threshold resonance this parameter can be smaller than that for the 
'?/'(3770) by a factor of 10 or more. In what follows it will be assumed that for the discussed 
resonance the value of w is between 1 and 10 MeV. 

As previously mentioned, for the actual D* mesons the isotopic mass difference gives rise 
to existence of two separate thresholds, the lower one for the neutral mesons and the upper 
one for the charged mesons. The energy difference^ 6 = 2M(D*+) — 2M(D*°) = 6.60 ±0.14 
MeV between these thresholds is generally not very small on the discussed energy scale 
and it is appropriate to take it into account. For the illustrative estimates in this paper 
we assume that the coupling of the resonance to the D* mesons is isotopically symmetric 
and also ignore the Coulomb effects, so that the only difference between the channels with 
neutral and with charged mesons is purely kinematical. Using the previous notation W for 
the excitation energy in the neutral channel: W = E — 2M{D*^), we write explicitly the 
production amplitudes in three separate regions of energy. Above both thresholds, i.e. at 
W>6: 

V2W-Wo + i [2ro + W ^ + {W - 6) ^/^] 

Affr = — ^ . , (6) 



V2W-Wo + i: 



A, 



ai + 



hL 



+ {W -5) 



W-& 



■ ^=T > (7) 

W-Wo + l [2V, + W + {W - 5) ^J^] 

where the subscripts n and c refer to the channels with the neutral and the charged mesons. 
Between the thresholds, i.e. at < < 5, the production amplitudes take the form 



Ahu — —:= 



hu k. 



H l^n 



US-W)J^+W-Wo + U2To + W 



(8) 



Aj 



ai + 



5-W 



+ W -W0 + U2T0 + W 



1(5 -W) 

Finally, below both thresholds, i.e. at W < 0, using the notation A 

bL 



(9) 



-W one can write 



A, 



{6 + A) 



S+A 



A 



(10) 
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Figure 2: The types of behavior of the e^e annihilation cross section (in arbitrary units) into 
one of the 'hght' channels near the D*D* threshold. The horizontal axis in each plot spans 
the cm. energy range from 3980 MeV to 4040 MeV. The D*oz)*o and D*+D*- thresholds 
are shown with shorter vertical tick marks. The plots are shown for three assumed values 
of the 'nominal' position Wq of the resonance X relative to the D*^D*° threshold: -20, 3, 
and 20 MeV, and the assumed value for this position is indicated in each plot by the longer 
vertical tick mark. In the plots in the upper row it is assumed that the coefficients ai and bi 
are relatively real and have the same sign. In the lower row of plots the relative sign of these 
coefficients is assumed to be negative. The solid lines in the plots correspond to w = 1 MeV, 
and the dashed one are for w = 10 MeV. The width parameter Fq is fixed at Fq = 50 MeV. 



That types of behavior of the cross section that can be observed in each of the 'light' 
channels, i.e. D*D, DD, or DgDs, is illustrated for some representative values of the pa- 
rameters. It is worth mentioning that the interference of the resonant amplitude with the 
non-resonant background is generally different for different channels, which may result in a 
significantly different behavior of the observed shape of the cross section for each channel. It 
is important to notice however that the features in the behavior of the cross section generally 
do not coincide with the position of the resonance. The shape of the cross section of the 
e"^e~ annihilation into the vector meson pairs D*D* above the threshold as found from the 
expressions (jHI) and (jHl) is only weakly sensitive to the parameters of the resonance and 
is dominated by the very rapid rise from the threshold due to the P wave factor k^. 



8 



Admittedly, at present the parameters involved in the discussed here description of the 
cross section are essentially unknown, and no prediction of specific shape of the cross section 
near the D*D* threshold can be made. However the presented here consideration illustrates 
that this energy region can be very feature-rich and a detailed scan of the e~^e~ annihilation 
cross section into each of discussed channels with charmed meson pairs can provide a very 
interesting information on the parameters intimately related to the strong interaction of 
heavy mesons, an area for which the information is quite scant. Even the determination 
of the basic parameters of the possible resonance considered here should necessarily involve 
the channel coupling effects. For instance, the width parameter Fq, taken in the plots in 
Fig. 2 at its PDG value [T] of 50 MeV can in fact be quite different due to the significant 
distortion of the cross section due to the discussed effects. An analysis of the angular 
distribution of the produced vector meson pairs would also determine the spin structure of 
the assumed threshold resonance. Finally, as mentioned previously, although the present 
discussion addresses mainly the D*D* threshold, similar phenomena may take place also at 
other heavy meson pair thresholds, e.g. near the threshold region for D*D pairs at bout 3875 
MeV, so that an exploratory test of those regions may also bring insight into new features 
of the heavy meson dynamics. 

I thank Dan Cronin-Hennessy, Brian Lang and Ron Poling for enlightening discussions 
of the CLEO-c experiment. This work is supported in part by the DOE grant DE-FG02- 
94ER40823. 
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